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Abstract 

The entanglement of graph states up to eight qubits is 
calculated in the regime of iteration calculation. The en- 
tanglement measures could be the relative entropy of en- 
tanglement, the logarithmic robustness or the geometric 
measure. All 146 local inequivalent graphs are classified as 
two categories: graphs vifith identical upper LOCC entan- 
glement bound and lower bipartite entanglement bound, 
graphs with unequal bounds. The late may displays non- 
integer entanglement. The precision of iteration calcula- 
tion of the entanglement is less than 10"^*. 

PACS number(s): 03.67.Mn, 03.65.Ud, 02. 10. Ox 
Keyword(s): graph state; iteration; multipartite entan- 
glement 

1 Introduction 

Multipartite entanglement plays an important role in 
quantum error correction and quantum computation. The 
quantification of multipartite entanglement is still open 
even for a pure multipartite state. Until now, a variety of 
different entanglement measures have been proposed for 
multipartite setting. Among them are the (Global) Ro- 
bustness of Entanglement [3] , the Relative Entropy of En- 
tanglement^ [2] , and the Geometric Measure 0]. The 
robustness measures the minimal noise (arbitrary state) 
that we need to added to make the state separable. The 
geometric measure is the distance of state to the closest 
product state in terms of the fidelity. The relative entropy 
of entanglement is a valid entanglement measure for mul- 
tipartite state, it is the relative entropy of the state under 
consideration to the closest fully separable state. 

The quantification of multipartite entanglement is usu- 
ally very difficult as most measures are defined as the solu- 
tions to difficult variational problems. Even for pure mul- 
tipartite state, the entanglement can only be obtained for 
some special scenario. Fortunately, due to the inequality 
on the logarithmic robustness, relative entropy of entan- 
glement and geometric measure of entanglement [5] [6] [7], 
these entanglement measures are all equal for stabilizer 
states [S] . It is known that the set of graph states is a sub- 
set of the set of stabilizer states. Thus these entanglement 
measures are all equal for graph states. 

The paper is organized as follows: In section 2, we 
overview the concept of graph state and graph state basis, 
the entanglement inequalities, the bounds of the entangle- 
ment. In section 3, we describe the iteration method of 



finding the closest product state. In section 4, we classify 
all the inequivalent graphs according to the entanglement 
bounds up to 8 qubits and propose the precise values of 
entanglement. Section 5 is devoted to the precision of the 
iteration. Conclusions are given in section 6. 

2 Preliminary 

2.1 Graph state 

A graph G = {V; T) is composed of a set ^ of n vertices 
and a set of edges specified by the adjacency matrix F, 
which is an n X symmetric matrix with vanishing diagonal 
entries and Tab = 1 if vertices a, b are connected and Tab 
— otherwise. The neighborhood of a vertex a is denoted 
hy Na = {v £ V \Tav = 1}, i-e, the set of all the vertices 
that are connected to a. Graph states [11] [12]] are useful 
multipartite entangled states that are essential resources 
for the one-way computing [13] and can be experimentally 
demonstrated [H]. To associate the graph state to the 
underlying graph, we assign each vertex with a qubit, each 
edge represents the interaction between the corresponding 
two qubits. More physically, the interaction may be Ising 
interaction of spin qubits. Let us denote the Pauli matrices 
at the qubit a by Xa ,Ya, Za and identity by la ■ The graph 
state related to graph G is defined as 



\G)= II Uab 
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^(-l)^^r^"|M), (1) 



where \fi)^ is the joint eigenstate of Pauli operators Za 
{a e V) with eigenvalues (—1)''", |+)^ is the joint -1-1 
eigenstate of Pauli operators Xa { a £ V) , and Uab {Uab = 
diag{l, 1, 1, —1} in the Z basis) is the controlled phase gate 
between qubits a and b. Graph state can also be viewed 
as the result of successively performing 2-qubit Control-Z 
operations Uab to the initially unconnected n qubit state 
\-\-)^ ■ It can be shown that graph state is the joint -1-1 
eigenstate of the n vertices stabilizers 



Ka^Xall Zb -.^XaZM^, O £ V. 



(2) 



beNa 



Meanwhile, the graph state basis are \Gki,k2.---kn 
riaey^a'-IG), with fc, = 0,1. Thus 

Ka \Gki,k2,---k„) = (-l)'"" |Gfci,fc2, •••*:„) • 



(3) 



Since all of the graph basis states are local unitary equiv- 
alent, they all have equal entanglement, so we only need 
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to determine the entanglement of graph state \G). Once 
the entanglement of a graph state is obtained, the entan- 
glement of all the graph basis states are obtained. 

A widely used local operation in dealing with graph 
states is the so-called local complementation (LC)|15|. 
which is a multi-local unitary operation Va on the a — 
th qubit and its neighbors, defined as Va = \/Ka = 
exp(— i jXq) JljgjY exp(i|Zfc). LC centered on a qubit 
a is visualized readily as a transformation of the subgraph 
of a-th qubit's neighbours, such that an edge between two 
neighbours of a is deleted if the two neighbours are them- 
selves connected, or an edge is added otherwise. 

2.2 Entanglement inequalities 

The global robustness of entanglement [3] is defined as 

R{p) = mint (4) 

such that there exists a state uj such that (p -|- tLo)/{'^ + t) 
is separable. The logarithmic robustness [9j is 

LR{p) = \og^{l + R{p)). (5) 

The relative entropy of entanglement is defined as the " dis- 
tance" to the closest separable state with respect to the 
relative entropy [2], 

Er{p) = mill S {p\\ uj) , (6) 

where S {p\\ uj) = —S[p)— tr{plog20j} is the relative en- 
tropy, S{p) is the von Neumann entropy, and Sep is the 
set of separable states. The minimum is taken over all 
fully separable mixed states uj. 

The geometric measure of entanglement for pure state 
IV'), is defined as 

EM)^ jmn -log^K^IV)^ (7) 

I (p) tz Pro 

where Pro is the set of product states. 

Hayashi et al [S] has been shown that the maximal num- 
ber N of pure states in the set \i = l,...,iV}, that 
can be discriminated perfectly by LOCC is bounded by 
the amount of entanglement they contain: 

log2 N < n-LR{\i',)) < n^EAm) < n~Eg{\^,)), (8) 

where n = \og2 Dh, Dh is the total dimension of the 
Hilbert space, and x = denotes the "average". 

For stabilizer state \S), it has been shown that [5] 

LRQS)) ^ Er{\S)) = E,{\S)) 

The technique is to prove LR{\S)) ^ Eg{\S)) and utilizing 
inequality ©.The entanglement can be written as 

Ei\S))^min-log,\{S\(b)\\ (9) 

where |(/)) = l\ji\/Pj\0) + \/^ - Pj^^"^' |1)) is the product 
pure state. 



2.3 Entanglement bounds 

The entanglement is upper bounded by the local operation 
and classical communication (LOCC) bound Elqcc = 
n — log2 N , and lower bounded by some bipartite entan- 
glement deduced from the state, that is, the 'matching' 
bound Ebi [TO]. It is well known that all graph states are 
stabilizer states, so the inequality for the entanglement of 
a graph state is 

Eu<E< Elocc- (10) 

If the lower bound coincides with the upper bound, the 
entanglement of the graph state can be obtained. This is 
the case for '2-colorable' graph states such as multipartite 
GHZ states, Steane code, cluster state, and state of ring 
graph with even vertices. For a state of ring graph with 
odd n vertices, we have [^\ < E < [f ] [TOj. 

The fidelity = \{G\ (^)|^plays a crucial rule in calcu- 
lating the entanglement. For a graph state, we have 

E = min -log2 |(G| cjj)\' = -log2( max F^). (11) 

Denote F = max^gp™ -F^ as the fidelity between the graph 
state and the closest pure separable state. The upper 
LOCC bound for a graph state is 

E<n-\A\, (12) 

since the largest number of entanglement basis states is 
2l^l, where \A\ is the largest number of vertices with any 
two of the vertices being not adjacent [TO] . 

The entanglement is lower bounded by the entanglement 
of a bipartition of the graph. The lower bound can be 
found by " matching" [TO] . A convenient way of finding the 
lower bound of the entanglement is to find the largest set 
of non-adjacent edges first, then assign the two vertices 
of each edge to two parties to form a bipartition of the 
graph. The vertices that are not assigned can be assigned 
to either parties. To verify if these edges are the last Bell 
pairs that can be obtained, one can apply local Control-Z 
and LC to delete the redundant adjacent edges. It can be 
verified that all the graphs up to 8 qubits in the literatures 
can be treated in this manner to obtain the lower bound 
of the entanglement. This is not a difficult task since most 
of the graphs in the literatures are already in the simplest 
LC equivalent form. Thus, at least for graph states up 
to 8 qubits, the lower bound of the entanglement can be 
obtained by counting the largest number of non-adjacent 
edges. 

3 Iterative method for the closest 
product states 

If the upper LOCC bound coincides with the lower bipar- 
tition bound, the entanglement of the graph state can be 
determined and equals to the bounds. Still there are graph 
states that the two bounds do not meet. We need a sys- 
tematical method to calculate the entanglement of such 
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graph states according to Eq. pT|) . The product pure 
state 1^) can be denoted as 

\(l>)^l[{x,\0)+y,\l)) (13) 

j 

where Xj and yj are complex numbers subjected to \xj\^ + 
\yjf = 1. Denote / = (G| <j)) , then 



(14) 



Let L — l/l^ — Aj(|xj|'^ + \yjf — 1), where Xj are 



the 



Lagrange multiphers. Then — and — lead to 
df 



dV3 



r - 



0, 
0. 



The two equations are combined to 
df , df 



y 



■' dxj ■' dyj 



(15) 
(16) 

(17) 



The left hand of Eq. ^ is (G| 0j) , with 



= l[(^k\0)+ykmy*\0)-x*\l)) 

n 

X Y[ (a;™|0) +2/™|l)). 



(18) 



Thus Eq. (4) is to say that the graph state is orthogonal to 
all {j = 1, . . . ,n) when |0) is the closest product state. 
It is clear that \<j)) is orthogonal to all {j — 1, . . . ,n ) 
too. Denote Zj — yj/xj, the derivatives are 



dl 



1 ^(_l)Kr.-J|(4-..,..) 



(19) 



n 



(20) 



where i/ = {fii, . . . , 0, ^j+i, . . . ,Mn}, 

i/'= {^1, . . . 1,/ij+i, . . . ,^„}, and 1' = 

{1,...,1,0,1,...,1}, 0' = {0,. .. ,0,1,0,. ..,0}. The 
binary vector 1' has all its entries being 1 except the 
j — th entry being 0. 0' is the logical NOT of 1'. From 
Eq. ()17p we obtain the iterative equations for Zj, 



k^j ^k 



k^J 



(21) 



We consider the change of fidelity in one step of iteration, 
that is, we only renew zj according to Eq. (PT|) while 
keeping all the other Zk {k ^ j) invariant in the step. Let 
^3 ~ Sx'' 9] ~ then hj and Qj are invariant in the 
step, / = Xjhj + yjgj- Forget the iteration equation for a 
while, we seek the maximization of the fidelity with respect 
to Xj = cos 9, yj — sinOe'^''^. The maximal fidelity should 
be \f\' 



\hj\^ + \gj\^ , which is achieved when 



El 



(22) 



The condition ([22]) of maximal fidelity is just the iterative 
equation pTjl . Thus in each step of the iteration, the fi- 
delity does not decrease. The fidelity increases abruptly 
or keeps unchanged in one step. In fact in each step, the 
fidelity can increase continuously from its initial value to 
its final value by changing {9, ip) continuously. 

Starting with any initial complex random vector 
z =(zi, . . . , z„), the iterative equation renews each Zj suc- 
cessively, the fidelity increases (or does not change) . After 
all Zj are renewed, a new round of iteration starts. The 
whole picture of iteration can be seen as a discrete process 
of the fidelity. The fidelity increases in each step until it 
does not increase any more. There may be the case that 
the fidelity reaches its local maximum. To find the global 
maximum, we run the iterative algorithm many times with 
random initial z. Moreover, we calculate the fidelity of the 
graph state with respect to random separate states for a 
million times to determine roughly the possible range of 
the fidelity before the iteration calculation. 

4 Classification of the graph state 
up to 8 qubits 

The LC inequivalent graphs up to 7 qubits are all plotted 
in |11| and numbered. There are a two qubit graph (No.l), 
a three qubit graph (No. 2), 2 of four qubit graphs (No. 3 
and No.4), 4 of five qubit graphs (No. 5 to No. 8), 11 of 
six qubit graphs (No. 9 to No. 19), 26 of seven qubit graphs 
(No.20 to No.45). In [H], the authors plotted all 101 LC 
inequivalent graphs of 8 qubits. The graphs of 8 qubits 
are numbered from No. 46 to No. 146. 

4.1 Graph states with equal lower and up- 
per bounds 

The entanglement of graph states with equal lower and 
upper bounds can be calculated with the methods in Ref. 
[TU] . It is listed in Table 1 and Table 2 for completeness. 

The graphs that are " 2-colorable" up to 8 qubits are 
listed in Table 1. It is a well known fact in graph theory 
that a graph is 2-colorable iff it does not contain any cycles 
of odd length. The LOCC upper bound and the lower bi- 
partite bound of the entanglement for " 2-colorable" graph 
state can be obtained by the methods described in Ref. 
[TU] . For each of these "2-colorable" graph states, it has 
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Table 1 



E 


No. 


1 


1,2,3,5,9,20,46 


2 


4,6,7,10-12,15,21-24,31,47-51,69-70 


3 


13-14,18,25-30,34,38,43,52-63,74,76-77,81, 


3 


83-84,103-104,122 


4 


64-68,87,89,91,95,99-100,120,128,143 



Table 2 

E II No. 

~3~' 16-17,32-33,35-37,71-73,75,78-80,82,102,121 
4 86,88,90,92-94,96-98,106-119,123-127,129-132, 
4 II 135, 136*, 144 

*the LC equivalent of No. 136 



been found that the two bounds coincide with each other 
Elocc — Ehi — Er, and the relative entropy of entan- 
glement is equal to the entanglement in Schmidt measure 
Er = Es = E. Table 1 shows the results. 
The LOCC bound for a "non 2-colorable" graph can be 
obtained with the largest set of non-adjacent vertices. The 
lower bipartition bound can be found by first searching for 
the largest set of non-adjacent edges, then verifying the 
candidate Bell pairs with local Control-Z and LC. Thus 
for "non 2-colorable" graph, the entanglement bounds of 
graph state can be obtained with "balls" (vertices) and 
"sticks" (edges) in the graph. When Elocc — Eu = Er, 
the graph states are shown in Table 2 (with E = Er = Es). 
No. 101 graph is special for the graph state has the relative 
entropy of entanglement Er — 4, the Schmidt measure 
Es = 3. 

4.2 Graph states with unequal bounds 

Up to 8 qubits, what left are the "non 2-colorable" graph 
states whose upper entanglement bound Elocc {Eu) and 
lower bound E\ji (Ei) do not coincide. Wc utilize Eq. ([2T|) 
to iteratively calculate the entanglement and find the clos- 
est product state with random initial complex numbers for 
Zj (j = 1, . . . , n). The values of relative entropy of entan- 
glement are listed in Table 3. 

A detail comparison of computed closest product 
states of No.8, No.39, No.41, No.45, No.85, No.l05, 
No.l34,No.l37,No.l38, No.l40 shows that aU these clos- 
est states have a substructure of the closest product state 
of ring 5 graph (No.8), although graph No. 140 docs not 
contain ring 5 graph explicitly, graph No.45 seems to con- 
tain graph No. 19 as its subgraph. Ring 5 graph is essential 
to all these graph states with entanglement A: -I- 0.9275 (in- 
teger k). In Ref. [TT] an identical product closest state is 
supposed for ring 5 graph state, and it has been shown 
that the entanglement of ring 5 graph state is 

Er^na5 = 1 + logs 3 + log2(3 - %/3) w 2.9275. (23) 
Denote |$,) = ^\0) + VT^e^'^^ |1) ,(j = 1, . . . , 4), with 
^ = ^5(1 -75) ~ 0.4597, ^1 = f,^2 = -f,(P3 = 



^ , (/34 = ~ X ■ Typically, the closest product state of ring 
5 graph state is 

\<f>r^ns5)^\'^lf\ (24) 

The other graph states may have their closest product 
states 

I0JVO.39) = |-)|0)|<I>3)|<I>2)^', (25) 

l^NoAl) = |<i>4)|<fl)^'|<l>4)|0)|-), (26) 

|</'^o.45) - |-))|0)|$4)|$3>'''|<I'4>, (27) 

|</>^o.85) = hf'|$l)^'|$4)|0), (28) 

\<Pno.105) = |+)^'|$l)^'|l), (29) 

|0iVo.l34) = |-)^'|<I'4>|$3)^'|$4)|0), (30) 

|0iVo.l37) = |+>|0>|$4)|$1>^'|$4>|0), (31) 

\cbNo.l3s) = |<1>2)|$3)'''|0)|+)|0)|$3)''', (32) 

|0iVo.l4o) = |l)|<f2)^'|<l>4)^'|0)h), (33) 

where |±) = ^(|0)±|1)). 

The next graph set (No. 19, No. 139, No. 141 ) with non- 
integer entanglement ( fc -I- 0.5850) graph states is specified 
by No. 19 ([[6,0,4]] stabilizer state). The entanglement of 
No. 19 is 

= 2 + log2 3 w 3.5850. (34) 
Typically, the closest product state is 

I'/'^o.ig) = |<i>3)'''|$4)''', (35) 

The closest states for No. 139 and No. 141 graph state can 
be 

10^0.139) = |-)|1)|<I'4)|<I>3)®'|$4)|<1'1). (36) 
\<t>No.l4l) = |-)|0)|<i>2)|<i>4)^'|$3)|<f2). (37) 

No. 139 and No. 141 graphs have No. 19 as their subgraph. 
The entanglement of No. 133 graph state is 

Eno.133 = 2 + 31og2 3 + log2(2- VS) (38) 
« 4.8549, (39) 

its closest product state can be 

|0iVo.l33) = |'i>4) \^lf^ l^if |$4) . (40) 

The entanglement of No. 44 is i?jvo. 44 — 4, its closest 
product state can be 

\cl^NoA,) = \Of'\ + f'\0)\-). (41) 

where |0) = ^(l*-*) ~* The closest product states for 
No.40,42,142,145,146 can also be obtained, the iteration 
calculation should be modified as explained in the next 
section. 
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Table 3 



No. 


Eu 


El 




Es 


Ps 


8 


3 


2 


2.9275 


2-3 


0.997 


19 


4 


3 


3.5850 


3-4 


1.000 


39 


4 


3 


3.9275 


3-4 


0.790 


40 


4 


3 


4 


3-4 


0.241(1) 


41 


4 


3 


3.9275 


3-4 


0.264 


42 


4 


3 


4 


3-4 


0.950(1) 


44 


4 


3 


4 


3-4 


0.432 


45 


4 


3 


3.9275 


3-4 


0.967 


85 


4 


3 


3.9275 


3-4 


0.689 


105 


4 


3 


3.9275 


3-4 


0.658 


133 


5 


4 


4.8549 


4-5 


0.969 


134 


4 


3 


3.9275 


3 


0.646 


137 


5 


4 


4.9275 


4-5 


0.917 


138 


5 


4 


4.9275 


4-5 


0.617 


139 


5 


4 


4.5850 


4-5 


0.999 


140 


5 


4 


4.9275 


4-5 


0.571 


141 


5 


4 


4.5850 


4 


0.935 


142 


5 


4 


5 


4-5 


0.281(1) 


145 


5 


4 


5 


4-5 


0.870(2) 


146 


5 


4 


5 


4-5 


0.501(1) 



Eu = Elocc, El = Eti 



5 Precision of iteration 

We concentrate on the precision of iteration for calculat- 
ing the entanglement of graph state whose lower and upper 



bounds do not meet. Let A = \E, 



numeric 



Ethe 



be the 



computational error of the iteration, where Enumeric is the 
entanglement determined by iteration, Etheory (= Er) is 
the entanglement proposed in the former section . We use 
the exact value of Etheory rather than its approximation. 
For simplicity, we just give the successful probabilities of 
achieving the precision within A < 10^^^ for some reason- 
able rounds of iteration with random initial conditions. 
From the actual numerical calculations, we can see that 
a precision of 10"^** is limited by the computer for our 
iterative algorithm (without double precision calculation) . 

For all graph states presented in Table 3 except No. 40, 
42, 142, 145, 146, the algorithm can be applied directly. 
The successful probabilities {Ps ) are listed in Table 3. 
The round of the iteration is set to 150 except for No. 140, 
whose round of iteration is 300. We renew z after each 
round instead of renewing Zj after each step of the round 
in the actual calculation for the reason of programming. 
In order to calculate the successful probability, we run the 
algorithm 1000 times for each graph state to count the 
number of algorithm that achieves the precision within 
A < 10-". 

For No. 40, 42, 142, 145, 146 graph states, direct appli- 
cation of iterative algorithm fails. The numerical results of 
entanglement are all greater than the values given in Table 
3, but the precision of the calculation is far from satisfac- 
tory. The precision can be 10"^ or so. A detail analysis of 
the separable state which gives best numerical value of en- 



tanglement shows us that the iterative equations ipij) are 
correlated. The common figure of these nonlinear correla- 
tions of equations can be illustrated by applying iterative 
algorithm to the simplest graph state, the No.l graph state 
(Bell pair). The iterative equations should be 



1-Z2 



Z2 



1 +21 



(42) 
(43) 



Substituting Eg. ((43)) into Eq. (|42)) . we obtain the identity 
z\ ^ z\. Thus the two equations are correlated. The cor- 
relation of equations leads to the fail of iteration. We can 
delete one of Eq.(f42|l into Eq. (|43|) to solve the problem. 
The fidelity is 



i/r 



II 



Zl + Z2 - 21^2 



4(l + |zi|^)(l + |zi|^) 



(44) 



Applying Eq. ([42|) and ignoring Eq. ((43|) , we obtain the cor- 
rect maximal fidelity |/| = \- Thus to obtain the maxi- 
mal fidelity, we should omit some of the equations and use 
the remain equations for iteration. For No. 40, 42,142 and 
146, we omit one of the equations, indicated in Table 3 
with notation (1) behind the successful probabilities. For 
No. 145, we omit two of the equations, indicated in Table 
3 with notation (2) behind the successful probability. In 
the numerical calculations, we set one or two zi to random 
numbers that do not change in the iteration, respectively. 
Since we do not know if all the equations are correlated 
or only some of them are correlated, we calculate all pos- 
sible choices of fixing Zi. For a given graph state, some 
of the choices of fixing zi may not lead to sufficiently high 
successful probabilities or simply fail. The successful prob- 
abilities shown in Table 3 are the best. 

We can see that the entanglement of all graph states in 
Table 3 can be efficiently calculated by iterative algorithm 
with very high precision. Most of them can be calculated 
directly, five of them can be calculated with modified iter- 
ative algorithm. 

A heuristic point of view is that we can set the fixed Zi 
to be 0. For an n vertices graph state |G„), the closest 
separable state should be |(/>„) = |(/>„_i) |0) when we set 
z„ = without loss of generality. Denote the n bit binary 
vector [In as (/i„_i, /^„), and the nx n adjacent matrix r„ 

r„-i 

c 



,then 



i(^„_i,0)r„(^„_i,0)^ = ^Mn-ir„-iM^_i. (45) 

Since ((/)„ |/i„) = {(t)n-i \ (0| /i„) = {(t>n-i \ (^o/^^, from 

the definition of graph state, we have 

(G„|<^„_i>|O>--^(G„„i|0„_i). (46) 

Where Gn~i is the subgraph of G„. From Eq. fi^ . a gen- 
eral relation for entanglement of any graph state and its 
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subgraph state follows 

En < + 1. (47) 

The equality holds for the case when = I'/'n-i) |0) 
is the closest separable state. For graph states No. 40, 
42,142,145,146, we can choose \4>n) = \4>n~i) |0) as the 
closest separable state when |(/)n_i) is the closest sepa- 
rable state of the subgraph state, the calculation of the 
entanglement can be reduced, we have 

En - En-i + 1. (48) 

for these graph states. The subgraphs of No. 40 and 
No.42 belong to the set {7Vo.l3, No.lA, No.l7, 7Vo.l8},the 
entanglement of all the subgraph states is 3. Thus 
the entanglement of No. 40 and No.42 is 4. The sub- 
graphs of No. 142, No. 145 and No. 146 belong to the set 
{iVo.40, 7Vo.42, 7Vo.44},the entanglement of all the sub- 
graph states is 4. Thus the entanglement of No. 142, 
No. 145 and No. 146 is 5. 

6 Conclusions 

The entanglement of graph states measured in terms of the 
relative entropy of entanglement (also with logarithmic ro- 
bust, the geometric measure) is given up to eight qubits. 
We use the iterative method to calculate the fidelity of 
graph state with respect to closest separable pure state. 
The iterative equations are results of the maximization 
of the fidelity. The equations have a clear meaning that 
graph state should be orthogonal to all other separable 
states that are orthogonal to the closest separable state. 
We have proved that in each step of a round of iteration 
the fidelity does not decrease. The iteration calculation 
has a very high efficiency if the resultant closest state does 
not contain |0) or |1) at all qubits (No.8, No.l9, No.l33). 
The precision of the iteration calculation can be less than 
10~^^ for all graph states up to 8 qubits with unequal 
lower and upper bounds of entanglement. To avoid pos- 
sible missing of the global maximum, with random initial 
parameters, we calculate the fidelity for each graph state 
a million times without iteration to determine its rough 
range, and calculate the maximal fidelity 1000 times with 
iteration. Iterative method brings us with the exact en- 
tanglement value of the graph state if we substitute the 
numerical closest separable state with its nearest exact 
one. For a given graph state, there are many local equiva- 
lent closet separable states, they all lead to the same exact 
value of the entanglement. The precision of the numeri- 
cal calculation is defined as the difference of the numerical 
and the exact entanglement. For some of the graph states, 
the iterative equations may correlate with each other. We 
analyze the situations and present a revised iterative al- 
gorithm to obtain the entanglement. In all the cases of 
unequal bounds, the entanglement may be equal to its up- 
per bound (integer) or in between the bounds (not to be an 
integer). For all non- integer entanglement cases discussed, 



we have found that the qubit states are the indis- 
pensable ingredients of the closest separable states. Based 
on our calculation, the non-integer entanglement graphs 
could be further classified according to the number of 
in the closest separable state. 
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